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ABSTRACT 



We consider N — 1 two-dimensional (2d) dilatonic supergravity (SG), 2d 
dilatonic SG obtained by dimensional reduction from N = 1 four- dimensional 
(4d) SG, = 2 2d dilatonic SG and string-inspired 4d dilatonic SG. For all 
the theories, the corresponding action on a bosonic background is constructed 
and the interaction with (dilatonic) Wess-Zumino (WZ) multiplets is pre- 
sented. Working in the lenge-N approximation, it is enough to consider 
the trace anomaly induced effective action due to dilaton-coupled conformal 
matter as a quantum correction (for 2d models s-waves approximation is ad- 
ditionally used). The equations of motion for all such models with quantum 
corrections are written in a form convenient for numerical analysis. Their 
solutions are numerically investigated for 2d and 4d Priedmann-Robertson- 
Walker (FRW) or 4d Kantowski-Sacks Universes with a time-dependent dila- 
ton via exponential dilaton coupling. The evolution of the corresponding 
quantum cosmological models is given for different choices of initial condi- 
tions and theory parameters. In most cases we find quantum singular Uni- 
verses. Nevertheless, there are examples of Universe non-singular at early 
times. Hence, it looks unlikely that quantum matter back reaction on dila- 
tonic background (at least in large N approximation) may really help to solve 
the singularity problem. 

PACS: 04.60.-m, 04.70.Dy, ll.25.-w 
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1 Introduction 



It is widely known that models of two-dimensional (2d) dilaton gravity 
may qualitatively describe the main properties of gravitational collapse and 
the evolution of the Universe taking into account quantum matter back re- 
action. This is due to the fact that reduction of 4d gravitational theories 
with matter often leads to effective two-dimensional (dilatonic) gravity with 
dilaton coupled matter. In a recent work [|l|], we studied the quantum cos- 
mology of four- dimensional (4d) reduced Einstein gravity with minimal 
scalars (working in large and s- waves approximations). As a result of the 
reduction, one finds 2d dilatonic gravity with dilaton coupled scalars and 
quantum cosmology (of 4d Kantowski-Sacks form) where dilatonic coupling 
may be interpreted as a second scale factor of the singular, non-singular or 
big crunch type |^. 

It is quite interesting to generalize such a study for supergravities (SGs) 
which will be the purpose of this paper. First, one can understand which 
different quantum SG may lead to which early Universe model in compar- 
ison with purely bosonic versions. Second, some SGs represent low-energy 
effective actions of superstring theory and so such cosmologies (limiting to 
our approach) may possess some characteristic stringy cosmological features. 
Third, such study may help to select the "best" theories in the cosmological 
sense from among the existing variety of SGs. 

The present work is organized as follows. In the next section we review 
the construction of A^ = 1 2d dilatonic SG with A^ scalar supermulti- 
plets. Section three is devoted to A^ = 1 Callan-Giddings-Harvey-Strominger 
(CGHS) SG with dilaton coupled scalar matter superfields. Working in the 
large A^ approximation, we take into account quantum matter effects and 
write the corresponding equations of motion (with non-zero matter fermions) 
on a bosonic background. The numerical study of such equations show 2d 
cosmologies where the Universe expands first and then shrinks, or where 
the Universe oscillates. In section four, we apply the same technique to 
study 2d dilatonic models with supermatter which is the reduced version 
from 4d A^ = 1 SG with matter. Then two-dimensional cosmology with 
a dilaton may be interpreted as a 4d Kantowski-Sacks Universe. Numeri- 
cal estimations show that it is singular 4d Universe in the situations under 
discussion. In section five, we work with A^ = 2 2d dilatonic SG with 
dilaton-coupled matter. The numerical analysis of two-dimensional cosmolo- 
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gies with the back-reaction of quantum matter shows a quahtative difference 
from the N = 1 SG case due to the presence of vector field in the bosonic 
background. Section six is devoted to string-inspired 4d dilatonic SG (its 
bosonic version) with N dilaton coupled matter multiplcts. Working in the 
large approximation and using the 4d anomaly induced action, we ana- 
lytically and numerically investigate equations of motion for conformally fiat 
4d cosmologies. In the cases under consideration we have found singular or 
non-singular at early times or trivial cosmological solutions. 
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2 General Action of = 1 Dilatonic Super- 
gravity 

The supersymmetric extension of the CGHS model |^ was recently consid- 
ered in . This theory is a partial example of more general two-dimensional 
= 1 dilatonic supergravity 0. The quantum corrections in the large N 
approximation for such theories with matter have been discussed in 0. Us- 
ing such quantum corrections one can also study Hawking radiation in 2d 
dilaton supergravity. In this section, we briefly review the general action us- 
ing the component formulation of ref . In = 1 dilatonic supergravity in 
this paper all the scalar fields are real and all the spinor fields are Majorana 
spinors. 

The general = 1 action of 2d dilatonic supergravity is given in terms 
^The conventions and notations are given as follows 

• signature 

• gamma matrices 

a b i:ab , ■ ab 

7 7 = + ^e 75 , 

CTafc = ^[7a,7b] = '^''"■blS ■ 

• charge conjugation matrix C 

Here the index ^ means transverse. 

• Majorana spinor 

1/; = ^'' = CiF . 

• Levi-Civita tensor 

,12 1 ab bo 

e — ei2 — i , e — — e , ^ab — —^ba , 
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of general functions of the dilaton C(0), Z{(j)), f{(j)) and V^(0) as follows 

+i [$ ® $ ® TpizmUy - [zm ® $ ® Tp($)]in^ 

+ [^('^•)]inv • (1) 

Here $ = (0, x, F) is a dilaton multipletQ and Sj = (aj,^j,Gj) is a matter 
multiplet. These multiplets have the conformal weight A = 0. VI^ is the 
curvature multiplet which is given by 

W = (5, r/, -S' -^-R- \rY%u + \r%) . (2) 
Here R is the scalar curvature and 

V = -IsY^.+ '-e-'e''^^,^,. , (3) 

oOf, = -ie~^eaf,e^''dxe2 - ^^-^TstVa • (4) 

The curvature multiplet has the conformal weight A = 1. Tp{Z) is called 
the kinetic multiplet for the multiplet Z = [ip, (, H) (in terms of superfield 
operators this corresponds to (V^ + XW)Z) and when the multiplet Z has 
the conformal weight A = 0, Tp{Z) has the following form 

Tp{Z) = iH,pCO^) ■ (5) 

The kinetic multiplet Tp{Z) has conformal weight A = 1. The product of 
two multiplets Z^. = (^9^, (k, H^) {k = 1, 2) with the conformal weight A^ is 
defined by 

Zi ® ^2 = {(pi(p2, + '^2(1, fiH2 + ^2Hi - C1C2) • (6) 



^The multiplet containing C{4>), for example, is given by {C{(t)),C'{(j))x,C{4>)F 
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This is familiar in the language of superfields as the components of Z\Z2. 
The invariant Lagrangian for the multiplet Z is defined by 



mv 



(7) 



corresponding to the superaction S = J d?9E~^Z. 

We consider only bosonic background fields below since this will be suf- 
ficient for the study of the cosmological problems under consideration. On 
the bosonic background where the dilatino x ^i^nd the Rarita-Schwinger fields 
vanish, one can show that the gravity and dilaton part of the Lagrangian 
have the following forms: 

(s^ + \r) - C'{cp)Fs\ , 



[$®$®Tp(Z($))]inv 
[Z($)®<|.®Tp($)]iny 

[^("^)]inv 
For matter part we obtain 



~ e 
~ e 



02n(z(0)) + 2Z'(0)0F2 

e[V\cl^)F + SVm . 



(8) 



N 



1=1 



® Tp(/($))]iny - [/($) ® ® Tp(S,)] 



mv 



N 



i=l 



+ total divergence terms . 



(9) 



In eq.(^ and "~" means that we neglect the terms containing fermionic 
fields except the kinetic term of ^j. Here we have used the fact that 

^^75^ = (10) 

for the Majorana spinors ^j. Using equations of motion with respect to the 
auxilliary fields S, F, Gi, on the bosonic background one can show that 

C"(0)r (0) - 2V{<P)Z{<P) 



S 

F 
G, 



C"2(0)+4C(0)Z(0) 
C"(0)\/(0) + 2C(<^)r((/.) 



C"2(0)+4C(0)Z(0) 







(11) 
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Especially for the supersymmetric extension ||5[ of the CGHS type model 

C{<f)) = 2e-2'^ , Z{(f)) = 4e-2'^ , = Ae'^^ , (12) 

we find 

S = , F = -\ , Gi = . (13) 
We should note > in the supersymmetric theory or unitarity is broken. 



3 = 1 Supersymmetric CGHS Model with 
Dilaton Coupled Matter 

After integrating out the auxilliary fields as we have shown in the previous 
section, the classical action of the CGHS type dilaton supergravity model 
with dilaton coupled supermatter in a purely bosonic background is given by 



Sr. 



d?xe 



^ i=i 



2n 



{r + H 



(14) 



As we are going to work in large N approximation we may take into account 
only matter quantum effects. Then the bosonic part of the trace anomaly 
induced effective action ^ (for a pure dilaton coupled scalar, see also 0) 
together with classical action has the following form[] 



S 



Sc + W 



d xe 



-2(f) 



271 



R + Ad, 



1 fA^ 1 A^^ 

— —R-R 

27r 1 32 □ 4 ^ 



1 ^ 

- t; E (^^.ci^^''a, + ^^YD^^^ + 4A 
i=i 

-R+^(j)R 
□ 4 



(15) 



^ The inverse of □ 



'd^^/^g^^^d^ is defined by 



°x ( -{x,y) 



Here x and y are the cordinates of the space-time and the subscript x oi O means the 
derivative with respect to x. In the conformal gauge (|T^), the ambiguity of the boundary 
condition is absorbed into the function t^(x^) which appears in (pT[). 
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In the (super) conformal gauge, 

^±^ = -^6^', ^±± = (16) 

the equations of motion given by the variations of g±±, g^-, 4> and Oj have 
the following form: 



= r±± 

N 

'-±± 



1 ^ 



i=l 

+ f |(5±05±0)P+ (^±05^0)} 

{-2a±pa±0 + di<p} + t±(a;±) (17) 

= 

= e-2^ (2a+a_0 - 4d+(j)d_(j) - X^e^^) + tI_ 

N N N 

--d^d^p - -9+(/>9_0 + (18) 

= e-2'^(^-4a+a_0 + 4a+0a_0 + 2a+a_p 

+ {-^a+(pa_</.) - ^a_(pa+</.) - (19) 

= ^+{e-"^^.x^) + ^-{e-^'^^+x^) ■ (20) 

Here t^(a;^) is a function which is determined by the boundary condition 
and T^^ is the energy-momentum tensor of the matter fermion defined by 

= 2^"'^ E + ^^l>^^,i^ ' 9,u^n'd,^i) ■ (21) 

i=l 

If we use the ^-equation of motion, we find 

tL = tL{x^), tI_ = 0. (22) 
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Since we are presently considering the cosmological problem, we assume that 
all the fields depend only on time t and replace d± ^dt. Then Eg . (|T7|) tells 
and t£±. are constants: = and = ^T. Then Eqs.(p!7|), ([18|) , 
(|19D and pO| ) can be rewritten as 

1 ^ 

= e-"^Udtpdtcl)-2idt<f)f) + -e-'^J2idta.y + NT 

-- {-2aipai0 + dU] + iVto (23) 
= e-2^ (29^0 - 49^09*0 - 4AV'') 

-f^^V -f 5*05*0 + (24) 
= e-2^ (^-49,20 + 4(9*0)2 + 292p + 4AV^-^^(9,a,)'j 

^9,(p9,0) - (25) 



2 4 

= dt{^-^'^dta,) . (26) 

Eq.(l26|) can be integrated to be dtai = Aie'^'^. Here Ai is a constant of the 
integration. Note that T can be absorbed into the redefinition of to 

to + T^to. (27) 



In the following we use the redefinition of 
Solving Eqs.(^) and (p^, we obtain 

5f0 = (3-^pe2^)(9,0)2-(2 + ^e2^)9,p9,0 + ^e2n9,p)2 

_^e^<^A2 + 2A2e2^ - -e^% (28) 
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1 + je^') p] 



_4e2<^ ( 1 + ^e^"^) A' -8(1 + -e'^ ] to + ^X'e'" 



(29) 



Here 



N 



Substituting (pSf ) and ([29| ) into (pSD we obtain 



(30) 



Ar2 

-lee-^"^ + 2Ar + — e^"^ (p + 1) 



-1 



X |f e-^^ - (8p + 4)e-^^ - iV (l + + ^e^^ (p + 1) p| {d, 
+ |-|e-^* + ^ (P + ^) + (P + 1)} ^^P^*'^ 



+ |-8 + ye^^ + ^e^nP + l)U^ 



+|4e-2<^-iV(p + l)}A'e 



(31) 



Eq. (|3T|) determines to from the initial conditions for 0, dt4>, p and dtp. Note 
that attempt to study quantum cosmology in CGHS-like model (working 
in s-wave approximation and sigma-model description with ad hoc choice 
of dilatonic potential) has been done in ref.0. In our impression, above 
approach is better for treating the dynamical dilaton effects. Moreover, it 
can be applied to wider class of models. 

Equations (|28|), (pQ]) and (pT)) can be solved numerically for several pa- 
rameters with the initial condition = p= ^ = ^ = Oatt = O. Typical 
graphs are given in Figs.l - 4. We also calculated the two-dimensional scalar 
curvature R, which is given by 



R = 8e-'^Pd^d- 



P 



(32) 
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In Fig.l, both the conformal factor p and dilaton field increase monotoni- 
cally with time and there is a curvature singularity in a finite conformal time. 
We should note that there was not found a solution where p monotonically 
increases with time in the bosonic model If p diverges near the singular- 
ity, the singularity will appear in the infinite future of the cosmological time 
defined by 

di = e^dt . (33) 

Note gtt = 1 when we use the cosmological time t. It should be also noted that 
there will not appear any singularities in the two-dimensional scalar curvature 
in the solutions of Figs. 2, 3 and 4. In Fig. 2, increases monotonically in 
time but p increases first and decreases later, which means that the universe 
expands first and then later shrinks. In Fig.3, oscillates and p decreases 
with small oscillation which means that an oscillating shrinking universe is 
realized. The scalar oscillating curvature goes to zero. The behavior in Fig.3 
is very similar to that found in case of the bosonic model reduced from 4 
dimensional model with Einstein-Hilbert action if we replace the two- 
dimensional curvature by the four-dimensional one. 

In Fig. 4, increases monotonically in time and p increases first but de- 
creases later, i.e., the universe expands first and then shrinks. The scalar 
curvature goes to vanishing. The behavior of solutions when = 100 as in 
Fig. 4 is not sensitive to the value of A^, which will imply that the system 
would be governed by the naive large- structure found previously. In the 
large limit, Eqs. (|23|) , (|2^) , (p5D are rewritten as 

= W + 1 - (9,pf ) 




_l{_25ip9i0 + 520}+to (34) 

= -^d^p-^dt<l)dt<l> + \d^<j> (35) 

= -^e''^A' + [~d,{pd,^)~ld^p] . (36) 

Especially in case of = 0, we can delete from Eqs.(|3^, (pSf ) and (^) 
and we find 

(1 + p){dtpf - 8top = (constant) . (37) 
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This tells us that there is a singularity when p = —1. This singularity 
corresponds also to a singularity in the scalar curvature R in (B^) since 



R = Se-^f^d+d-p ~ 
Using a new variable and the parameter 

P={l + p)'^ , 



2e2 (-8to + 
1 + P 
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we can rewrite (|37|) as follows 



-(d,pf-9to{p 



E . 



(38) 



(39) 



(40) 



Therefore we can regard the system as that of a particle with unit mass in 

2 

the potential —9to{p3 — 1). When to > 0, the system becomes unstable and p 
and p increase monotonically. The to > case corresponds to an expanding 
universe. On the other hand, when to < 0, the universe shrinks even if it 
expands for some initial period and then later the universe "encounters" the 
curvature singularity at p = —1. 

Before closing this section, we make some remarks about the duality struc- 
ture of = 1 dilatonic supergravity on bosonic background. The duality 



might be useful when we consider the so-called "graceful exit" problem |T0 
for inflationary universe. 

We now consider the following simplified classical action: 



S= — 



2tt 



d Xy 



-9 



N 



-2<t> 



i? + 4(V0)' 



2a<j) 



1=1 



(41) 



Here we introduced a parameter a and we have put 
anomaly induced effective action [R] is given by 



0. Then the trace 



W 



-- I d^Xy 



-9 



N 



N 



N 



32n A 



R—R - —a^V^(f)Vx(j)—R + —a 



An 



(42) 



In the conformal gauge, the actions 

1 



and 



take the following forms: 



S + W 



d X 



-N\—^pd^d.p-c?pd, 



+ acpdj^d-p 



(43) 
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Here we treat S + W as a classical action and we solve the equations of 
motion. If we are interested in the case where all the fields only depend on 
time t, we can replace d± by ^dt- Then the action (^) becomes 



S + W 



/ d'^x 

271 



a 



a 



+N\--pdtp--p{dt<i>y + - 



(44) 



Furthermore, we replace the conformal time t by the cosmological time i 
defined in Eq.(|33|). If we use a new variable R instead of p 



(45) 



we obtain 



S + W = -— / d'^x 
In 



1 r 



;.2 



4 r 



16 8 



1 

16 



a 

-in 



a 

'T 



Inr + 



(46) 



Here = ^ and = The variables (j) and r express the classical duality 
(j) ^ r. We now consider further redefinitions of the variables. When 

= 8 

defining the following new variables 
U = 



(47) 



+ ATe^-^ (-lnr + — 
V4 32 



1 

32 



V 



132 8 3^^ 3 



we obtain 



S + W 



1 

2^ 



d^xUV + 0{N^) . 
14 



(48) 
(49) 



This tells that the system has a duality 

U (50) 

when is large but Nh is small. Note that the Planck constant h is implicitly 
contained in (|49|), i.e., fi was chosen to be a unity and 0{N'^) means 0{h^N'^). 
In case of the bosonic model, similar argument holds and we can find a 
"perturbative duality" in (|50| ) if we change the coefficients and parameters a 
little bit, e.g. = y instead of (P). 



4 4d Reduced Supergravity Model 

We now restrict our discussion to 4D action of Einstein gravity with a 
cosmological term and matter described by minimal scalars aj coupled to 
the metric with spherical symmetry: 

ds^ = g^^dx^'dx" + e-^-^cifi . (51) 

Here the two-dimensional metric and dilaton depend only on time and radius. 
The spherically reduced action reads 



Sred = J d'^Xy/^e" 



-2(j> 



1 1 ^ 

.{R+2{V<pr-2A + 2e"^>} + -J2{^X^ 



16nG ' ' ' 2 . , 

2=1 



(52) 

In the 4-dimensional = 1 supergravity model, the metric ( ^T] ) can be 
realized by choosing the vierbein fields as follows 



ej = e*, e2=sin^e^ = el = . (53) 



The expression (|53|) is unique up to local Lorentz transformation. The local 
supersymmetry transformation for the vierbein field with the parameter ( 
and ( is given by. 

Here ip^ is the Rarita-Schwinger field (gravitino) and we follow the standard 



notations of ref.|TT[] (see also [12]. If we require that the metric has the form 
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of (|5lD after the local supersymmetry transformation, i.e. 



^9te = Sgre = ^Qt^ 
5g^^ = sin 66 gee , 



Sgrw = Sge 



we obtain, up to local Lorentz transformation. 



sin6'?/;ei , 
- sin 9'ijje2 



,1 — 



(55) 



(56) 



i>l = sin Oipl , 
, V'J = - sin OipQ , 
-ipei , ipe = -ipe2 , 

i^l = -2e~'^4iJe2 , i^ti - ^\ = -2e-'f>elilJe2 



ei ■ 



(57) 



Eq. (|56|) shows that the local supersymmetry of the spherically reduced theory 
is = 1, which should be compared with the torus compactified case, where 
the supersymmetry becomes N = 2. Let the independent degrees of freedom 
of the Rarita-Schwinger fields be denoted by 



2^1 = ^rl + 

2^l=ija + ijl 
Xi = i^ei , X2 = 



2i)l = il)r2 + t 



2 

r 1 



■^92 , 



i^t2 + i^t 



(58) 
(59) 



then we can regard ilJt,r and x as the gravitino and the dilatino, respectively, 
in the spherically reduced theory. 

We concentrate on the spherically symmetric metrics of Kantowski-Sacks 
form [0 in ( ^2]) where g^u = a?{t)r]^y. Such a metric describes a Universe 



with a 5*^ X S*^ spatial geometry (for some of its properties, see JT^)- In the 
case of dilaton supergravity obtained by reducing of 4d supergravity with 
chiral scalar supermatter, the action in the bosonic background is given by 



S 



d xe 



27rl32 □ 4 



-20 

2. '^ + 2^" 



- 2A^ + 2e^^ 



□ 4 



(60) 
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Here, the last line term is the anomaly induced effective action due to matter. 
We obtain the following equations replacing (|^), (|29| ) and (|3lD 



-^e^'t'A' + X'e^P + e2*+2^ - je^% (61) 



- I -lee-^"^ + 2N+ — e^*^ (p + 1 



X 



|e-^ - (8p + 7)e-* - f (p + 1) + ^e^* (p + 1) p| (9., 



+ <!2e-^^-^-^e^nP + l)[(5*p)' 



P 

AT iV2 



/V2 

^ ) + —^"^ (p + 1) ^ 5*p5.- 



+ je^' + ^e^'ip + l] 



A' 



+{4e-20 _ N{p + l)}^e'^ - (Ae-'^f + ^) 



(63) 



Equations (|6lD , (|6^ ) and (^) can be also solved numerically for several 
parameters with the initial condition (f) = p= ^ = ^ = Oatt = 0. Typical 
graphs are given in Figs. 5 and 6. We calculated the four- dimensional scalar 
curvature corresponding to the metric (pT|) which is given by 



(64) 



In the 4d reduced model, the curvature singularity always seems to appear, 
at least in cases under discussion. In Fig. 5, both the dilaton field (j) and the 
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conformal factor p increase monotonically in time as in Fig.l in CGHS type 
model. If we regard the solution as describing a universe of the Kantowski- 
Sacks form with the topology S"^ x S^, the radius of S"^, which is given by 
e""^, decreases to zero but the radius of S^, which is given by e^'', increases. 
Note that the radius of corresponds to the radius of the universe when 
we regard this model as a two dimensional one. The four- dimensional scalar 
curvature increases first and goes to minus infinity. In Fig. 6, increases 
monotonically and p decreases at first but increases infinitely after that. This 
means S"^ shrinks but the factor shrinks at first and later expands. The 
four- dimensional scalar curvature goes to zero at first but increases to infinity 
later. The dilaton field always increases and runs away to the singularity, 
which means 5*2 in the universe of Kantowski-Sacks form shrinks to a point, 
which causes the curvature singularity. When we regard the model as a two- 
dimensional one, the radius of the universe goes to infinity at the final stage 
in the cases studied here. We see that at t = there is no singularity in most 
of our cases. However, at late times there is as a rule singularity. Note that 
the behaviour of our cosmologies at late times is not so important as other 
effects should define late time structure of the Universe. 



5 N = 2 Dilaton Supergravity Model 

In this section, we discuss N = 2 dilaton supergravity. Notations and con- 
ventions are mainly following to [Q. However, the possibility of considering 
a purely twisted matter potential term point out in |^ is not treated. The 
general classical action of [/^(l) = 2 dilaton supergravity with matter is 
given by^ 

S = — I (fx(fe£-^C{^)R+- I (fxd^e^E'^Z{^)Z{^) 
2tx J 2 J 

1 M _ 

+ J d^xd'^0£'^V{<^) + -J d'^xd^9E-^g{<^)Y,^^^i 

N 

+ / d^xd^dE-^f{<^) Xi^i + h.c. 



i=l 



i=l 



^We use = 1 or = 2 to denote type of SG but also we use N as number of matter 
multiplets which should not cause confusion in the context used. 
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d xe 



C'{(f))FB + C'{(I))FB 



1 

+2(C(0) + C{(t))){BB - 7^) + 2z(C(0) - 



1=1 

+d+{g{(j))ai)d_di + d+aid_{g{(j))ai)} 

- ^ ( - 1 ^^ 

+ /(0)) ^ - - {^^X^^+X^ + ^+X^^^X^)\ 

+fermionic terms (65) 

Here $ = (0, x,F), Sj = {ai,^i,Gi) and A' = {xiXi^ Hi) are the dilaton, 
chiral matter and twisted chiral matter multiplets, respectively, TZ is the 
scalar curvature and JF = OqAi — diAo is the field strength of the vector field 
Afj^. The action (|65D shows that the imaginary part of C(0) can be identified 
with the axion field. Integrating out A^, we find 

^(0) - ^(0) = c (constant) (66) 

If we consider the solution where c = 0, the dilaton field becomes real 

= 0. (67) 

If we consider the case of ^ M, 1, the bosonic part of the effective action 
would be given by 

Here 

and we used that the action of the fermion component of Xi in the bosonic 
background (|65D has the following form 

Sc = Wc^'a;e-i(/(0) + /(0))E (^+C-)C + (5=C+)C+ 

i=l 

+^(A^C-C-^=C+C+)} (70) 
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The action ( |68D is obtained by replacing N by 2N in the = 1 supergravity 
effective action and adding the second term coming from chiral anomaly. This 
is because the dilaton dependent function + /(0)) in front of the matter 
fermion kinetic term in Eq.(|70D can be absorbed into the redefinition of ( and 
the matter fermions do not give any contribution to the dilaton dependent 
terms of effective action. The contribution of the energy momentum tensor 
from the matter Q fermions is, as in iV = 1 case, absorbed into the redefinition 
of to- Therefore the behavior of the cosmological solution is almost identical 
with the = 1 case in the bosonic background with JF = 0. 

We now consider the case where ^ 0. For this case, we parametrize 
the vector field as follows; 

16 

A, = d,9 + -g,,e''Pd,9 . (71) 

Note that 6 expresses the gauge degree of freedom. Then we can rewrite the 
chiral anomaly term in Eq. (|68|) in a local form, 
N 

The ^-equation of motion obtained from S + W has the following form 

= □ {2z(C(0) - C(0)) + e} (73) 

that is, 

9 = -22(C(0) - C(0)) + constant . (74) 
Now we consider the case, 

C(0) = -^e-2^ Z(0)=4e-^ V(0) = Ae-2^ f{<P)=e-"^ (75) 

and the background where the chiral matter multiplet Sj vanishes. Then 
integrating the auxilliary fields B, B, F and F and using ([7^), we can rewrite 
the bosonic part of S* + as follows: 

S + W 

= — I dPxe-^ -e-2^7^ - , ^ |4e-^^9.(^9^(^ + 8^98^9] 
2nJ [ + 9^ ^ my- a. / 

^ ^^9U9 + 2e-^-Y.^X^^'x^ 



I d\^:F^T Sx^9u9 . (72) 



-i(e-4^ + ^2)e2^ + 16Ve-^^ + 9"^ N .^^ 
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Here 6 should be understood not to represent the degrees of freedom of the 
vector field but to be the imaginary part of e"^"^ using (|74D i.e. to be the 
axion field. Then the equations corresponding to (|2^) - ( PB] ) become (T is 
absorbed in the redefinition of to) 

= 4e-'-d,pd,^ - ^^J^^ ^ (4e-^^ {d,^f + {0,6^) - 1 {d,ef 

-I N j\j / 1 \ 

+ 2^-'^ E 9tX^^tX^ + J {db - (dtp?) + ^ (P + 2 ) (^*^)' 
N 



- {-2dtpdt^ + + 2iVto (77) 
-j^'p-y(5*^)' + y5*V (78) 



(e-4v + ^2)f^ (e-4<^ + ^2); 

.2^2g2p r_ (e-^'^ + g^)e^^ ^ ^ ^2)1 



|l(e-4^ + ^2)2^2^ _ 32(g-4^ + ^2)ig-4^| 

1 ^ r /V 1 

je-^-J^ftXiftx. + |-iVa,(M'P) - jd^p] (79) 



-2 



{e-^^ + e^)l \(e-4^ + ^2)i AT 

^ye ^ 2 ^e ^ 

(e-4^ + ^2)f ^ (e-4^ + ^2)f *^ 

-32A2e2''^(e-^^ + 9')'^ |_ (e^^^^^V^^ ^ ^g^^_4^ ^ ^2^1 1 ^gQ^ 

= ^t{e-^'^^tx^) = ^t{e~'''^tx^) • (81) 
Eq.(|8T|) can be also integrated to be dtXi = Aie'^'^ {dtXi = ^je^"^). Defining 
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= m^^=i^i^i^ Eqs.(|73), (0), (|7|) and (|0|) are rewritten as follows, 



-2 1 



AT 



dtpdt^p 



-2ip 



N 



+ 



r + 2 pe 

(e-4^ + ^2)i 2^ 

AT 



2ip 



2A2(e-^^ + ^2)e2p+2^ 



2 +_(£ll^±£!M^ + l6(e-4^ + ^2)i 

16 



(82) 



N \ ^ 
9(1 + -e^'^ j e-2'^9<p9t<^ 



A^ 
+ I2 + 



AT 



2^ 1 + 



A^ 



A^ 



+ + ^e^^) (dtp)' - 4e2^ f 1 + ^e^^) A 



,4(e- 



(e-4^ + ^2)i 2 



■pe 



2ip 



N 



4 

4^2g2p+2^/g-4^^^2 



1 + 



4 

A^ 
— ( 
4 

4 



+ 16(e-4v + ^2)^ 



,2</. 



(83) 



'1 r 



(e-4v + ^2)i + AT] 
2ee~^^ ,0 N2 

-32AV^(e-^^ + ^')^ 
4iVe-2'^ 



e 



2e-^^ 



dtifdtO 



;e-4^ + ^2) 
' (e-4^ + ^2)e2^ 
^ 4 
Ar2 



+ 16(e-^^ + ^2)i 



(84) 

-2n 



-Alp 



- + 16e~2^ 1 e^^ - AT^e^^p 
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+ ^2) 

4(e-4^ + 2^2)e-4^ ^ ^^-^^ ^ _ ^ 



(e-4^ + ^2)^ 



+ ^2 
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2A^ - 8e~^^ 1 e^'^ 



i (6-4^ + ^2)^ \^ 16 



P 



+2iV (^1 + ^e^^) p^dtpOt^ 

\2(e-4^ + ^2)i V 16 y 8 ^/^*^^ 

r + . . + 



+ ^e-2^ f 1 - — e^^^ - Ne'^ ( + -] p]idt9f 

-dt9dt(p 



(c-t'^+6'2)e2¥' 



2A2e2p+2^(e-4^ + e^) ( ^<cl^ -Np + 4e-2^ - A^) 



-4ly3 



2^ , + i^e«^,_8^_^^e*^)_^,e^A 



+ ^2)^ 2 " V 16 



Equations (^) , (|83D , (|J) and (|85D can be also solved numerically for several 
parameters with the initial condition (p = p = ^ = ^ = 0, 6 = 1 and 
^ = 1.1 at t = 0. Typical graphs are given in Figs. 7 - 10. In most of cases, 9 
increases linearly with respect to the conformal time t. As in the 4d reduced 
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model of Figs. 5 and 6, the singularity always seems to appear, at least in the 
cases under discussion. As a specific feature in the N = 2 model, there is a 
case (Fig. 10) where the dilaton field ip goes to minus infinity. Such a case 
has not been found in the = 1 model of Figs.l - 6 but was found in the 
bosonic model When < in the bosonic model, the dilaton field (p 
goes to minus infinity but note that or is always positive or zero in the 
supersymmetric model. If we regard this N = 2 model as a four dimensional 
model with the topology S"^ x of the Kantowski-Sacks form, the radius of 
S^, which would be given by e""^, goes to infinity. 

In Fig. 7, the dilaton field ip increases monotonically in time and goes to 
infinity, the conformal factor p decreases at first but increases after that and 
goes to infinity and 6 increases linearly in the conformal time t. In Fig. 8, if 
increases monotonically in time and goes to infinity, p increases at first but 
decreases after that and goes to minus infinity and 6 increases linearly in the 
conformal time t. Fig. 9 is the case that 6 increases monotonically in time but 
not linearly with respect to the conformal time t. In Fig. 9, both of and p 
decrease at first and increase after that and go to infinity. This case should 
also be compared with that of the A^ = 1 models, in most these the dilaton 
field (fi increases monotonically in time except in the case of Fig. 3. Even in 
Fig. 3, (f cannot be negative. In Fig. 10, the dilaton field ip oscilates a little and 
finally goes to minus infinity, the conformal factor p decreases monotonically 
and goes to minus infinity and 6 increases linearly in the conformal time 
t. In the cases of Figs. 7 and 9, the radius of the two dimensional universe, 
which is given by e^, goes to infinity, that is, the universe expands. On 
the other hand, in the cases of Figs. 8 and 10, the two dimensional universe 
shrinks to a point. It is, however, not clear what determines the fate of the 
two dimensional universe since it seems that there is no any clear correlation 
between the behavior of p and the "cosmological constant" A^ and/or the 
matter energy density A^. We again see that at t = there is no singularity in 
most of our cases. However, singularity appears in course of evolution. We 
again may argue that presented effects are getting negligible at late times 
where quantum corrections almost play no role. 
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6 String inspired model in 4 dimensions 



The low-energy effective action of tlie dilaton-gravity sector in the super- 
string would be given by 



S = J d^x^^e-^-^* {R - 6V^0V^0*) + 



(86) 



Here • • • denotes the terms containing moduli, gauge fields, which depend 
on the detail of the compactification, and fermionic fields. A special 4d 
dilatonic supergravity coming from string theories has been considered in 



(recently it has been independently rediscovered in |]T6[). This model 
specifically (and is expected to be most directly related to superstring and 
heterotic strings) uses a real linear multiplet instead of the more familiar 
chiral multiplet description that follows. We now consider the system where 



the massless matter fields couple with the dilaton gravity system (pq). 

We first consider the case that the matter does not couple with the dila- 
ton. In this case, the action of the matter field (WZ action) is given by 



+ h.c. 



(87) 



Here we used the notations of ||TT| . Then the bosonic part of the action is 
given by 

= -eAA*n - ed^Ad"^ A 
6 

+ ^{MA* - 3F*){M*A - 3F) - ^eAA*bab'' 
-'-{A*dmA-Ad^A*)b^ . 



(88) 



Note that the coefficient ^ in the first term, which assures the conformal 
invariance. The action ( ^8|) is, in fact, invariant under the following (local) 
conformal transfomation with the parameter a: 



9mn 

M 



A 
F 



ymn^ 

e^e ) 
Me-'' 
6™e-" 
Ae"" 
Fe-^-^ . 



(89) 
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The terms containing fermion fields are given by 



r L 

U = — -e[xo"™PmX + X^^^mX] + terms containing gravitino . (90) 



Here 



-mn 



a 



[a, 
(a, 



m ^ ^mnl ) X 



nJ 

O- UJmnl ) X 



1 



4 
1 



— m -n — n „m\ 

a a — a a 



— [a a — a (J 



Under the conformal transformation (^), ujmni transform as 
Then under the conformal transformation fBSl) and 



the action 



transform as 



C^ + 6Cf 
5Cf cx (xa'"x + x^'"x)5mcr. 



(91) 



(92) 



(93) 



(94) 



Since X'^""^ = ~4'^"'Xy vanishes and is invariant under the conformal 
transformation. 

The conformal invariance holds if we introduce dilaton (chiral) multiplet 
^ whose first component is invariant under the conformal transformation 
and the last component G is transformed as 



G 



Ge 



(95) 
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as follows 

£ = 1^2922^ --^(P:D-8i?) ($+(/(*)$) + (/(*+)$+)$) 
= ^e{f{<P) + f{<P*))AA*n 

-^e(/(0) + + Xct'^V^x] 

+—{{MA* - 3f{(l)*)F* - 3f'{(l)*)G*A*){M*A - 3F) 
18 

+ (MA* - 3F*){M*A - 3/(0)F - 3f\(p)GA)} 
-^e{f{<p) + f{<P*))AA*bab' 

+ ifi(l)*)A*dmA - Admifi(l>lA*))b"'} . (96) 

Here we put the fermionic component of ^ to vanish since we now regard ^ 
as a background. 

If we consider the background where the dilaton field is real, the non- 
local effective action W induced by the conformal anomaly of the dilaton 
coupled matter fields has been calculated in ref. []T7| and is given by: 



W = bj d^xy^Fa 
+b' J d^x^[ 



a 



a 



1^ ^" + + ^'0 / "^'""^ " " 6(Va)(Va)]' 
+ / rf^xv^{^[(V^)(V¥.)]V+|n((V¥.)(V^))a 

+|(V(^)(V(^)[(Va)(Va)]} (97) 

Here the a-independent terms are dropped, the last terms represent the con- 
tribution from the dilaton dependent terms in trace anomaly. The dilaton 
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dependent function (p is defined by 



(98) 



The coefficients are given by 
1 



120(47r) 



xAN , b' 



1 



ai 



32(47r)2 



x2N , 02 



360(47r)2 
1 



X 



24(47r)^ 



13A^ 

X 2N 



(99) 



Tfie coefficient 6" can be changed by a finite renormahzation of local coun- 
terterm and we can put b" = (as is obtained also from direct calculation). 
We consider below N WZ multiplcts and apply the large expansion which 
is why all the coefficients above have the multiplier N. In the following, we 
only consider the case where f{4>) is given by 



(100) 



and we assume (p is real and that the conformally fiat fiducial metric is 
Qfiu — G^'^Vniy- Then the total action has the following form: 

S = W + S,i 

= J d*x{2b'{aaf - [3b" + 2{b + b')] [aa + {d^af]" 

+|[a,^9V]V+|n(9,(^9V))<T 

+ ^df^(pd^(pdi,ad''a 

-Ge^^-^^d^ipd^ip + e""-^^ [-GOa - 6{d^a){d''a)]} . (101) 
The equations of the motion given by the variation of a and are given by 
= 46'nV 

-2[36" + 2(6 + b')] [□ [aa + {d^af] - 2d^ {a'^a [na + {d^af] }] 
+^[d,<pd^<pr + ^n{d,<pd^<p) 



02 
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_2e2-2v[-6n(T-6(9^a)(9'^(T)] . 



(102) 



(103) 



Since we are interested in cosmological problems, we only consider the uni- 
form solution, where all the fields depend only on (conformal) time t and we 
obtain 







-2[36" + 2(6 + 6')] 



8 \dt Adt^ 




d'^cr I da^ 
'd^^X'dt 



^d\da 
dt I dt 



d^a I da^ 



dip\ \ a2 d i da I di^\ 



_24g4<7~2^ 



dip^ 
Itt 



+ 2e2'"~2^ 



2 dt \ dt \ dt 



_6— - q(— 
df^ V dt 



^2 2.-2^ d . .da' 

-6 -— + 12— e^'^-^'P — 

dt^ dt \ dt 




= ^^{{^ 



+12^ fe— ^Ul2e— 



a2 d I dLpd^a\ a2 d ( d(f /' da' 
Ydi \dt~dlF) ^ Ydi I dt \dt 

2 



_2g2a-2^ 



c?2(T / (icr^ 
' ~d^ ~ \~di 



(104) 



(105) 



We can often drop the terms linear to a in (rUT) [T^. (Actually, if one 
considers the correspondent theory as IR sector of quantum gravity |jl9[ or 
IR sector of induced supergravity |2^ the omitting of linear term on sigma 
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corresponds to trivial redefinition of tlie correspondent source.) In tliis case, 
(|104|) and (|105|) can be reduced to 







46' 



■2[3b" + 2{b + b')] 

a2 d I da f d(p\ 
"2di \~dt \~dt I 



'd^ 


d'^a 

[dt^ + 


(daV 


2 ^ < 


1 da 


'd^a 


dt^ 


^dt J 


di ^ 


\~dt 


dt^ 



da 
dt 



,fa 
' dt^ 



-6- 



d 



da" 



+ 12— I e'"-'^ — 
dt 



dt^ dt 

2 



(106) 



a2 d I d I da\ 
'Ydi [dt \dt) 



+12^ /'e4.-2.M+i2e^.-2. 
dt \ dt J 



'dip' 



-2e 



2a-2ip 



' da^ 
Ht 



(107) 



In the large N limit, where the contribution from the classical part can be 
neglected in Eqs. (|102| ) and (|103|) , some analytical solutions can be found. 
When the terms linear to a in ( |101| ) are dropped, such analytical solutions 
were found in |21[]. On the other hand, when we include the a linear terms, 
we can construct only a few closed form solutions. Some of them are given 
by 



(1) (p = ipQ (constant), a = 

(2) (p = (fio (constant), a = aln (^-^^ , 



= 1 



2b' 



3b" + 2(6 + b'] 



(3) 



(J = 0, 



«2 



= 12— > 



> 1 



(108) 



In the case that the terms linear to a in ( |101D are dropped in the large 
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limit, ( |106| ) and ( |107| ) can be integrated to be 

d^a f da\^ a2 d j f dip 
T—rir — 2q 



dt^ \dt J A dtWdt 

2 

= e (109) 



dip ( da 
dt \ dt 



where c and e are integration constants and r = —(36" + 26), q = —[36" + 
26' + 26] . We can find some particular solutions of Eqs. ( |109| ) . 

1. For r = 0, we get 

a = a^t + const 

( et\ 

f = const exp ( — 1 (HO) 



where 



026^ /ae^ c 



That is a singular solution, it has an initial singularity at physical time 
equal to zero. 



2. For r 7^ 0, we find the general solution: 



ydy ^ f dy e f dy 



u{,y) ' J u{y) ' 2 J y'^uiy) 



where 



^2 



u\y) = \(qy'-'^ + 2cy + e\ (113) 



r \ y 
One particular solution in ( |113| ) is given by 



p = —-t , a = t , c = —2q — 026^ (114) 
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That is the same solution as in ( |110| ). Another particular solution is 
given by g = 0, c = 0, e = 1, r = 26', and 



a = a2\/\2b' 



1 . /(|26'|)-V2t + a2' 
- arcsm 

2 I a2 



1/2 



1 ( f (\2b'\)-^/H + a2 

H — sm < 2 arcsm 

4 I V a2 , 



1/2- 



;ii5) 



and similar solution for (f. That is again singular solution. In general, 
dilatonic coupling acts against such non-singular solutions. We did not 
find any non-singular solutions from the general solution (|112|) . 



We now consider solving ( |104| ) and ( |105| ) or ( |106| ) and ( |107| ) numericaly. 
Some examples of the solution by numerical calculations are given in Figs. 11 
- 14. The results in case that we include the terms in ( |101| ) linear to a ( |104| ) 



and (|105|) are almost identical with the case that we drop the terms linear to 
cr ( |106| ) and ( |107| ), as expected. For this reason, we only give the results in 
Figs. 11 - 14 for the case that we drop the terms linear to a. In Figs. 11 - 14, 
we calculate (solid line), a (broken line) and 4d scalar curvature (dotted 
line), which is given by 



R 



-6 




(116) 



for several and for the several initial conditions. An interesting thing 
is that the 4d scalar curvature oscillates as a trigonometric function. In 
Fig. 11, we choose A^ = 1 and the initial condition 0(0) = cr(0) = 0, 0(0) = 
a(0) = a(0) = (t(3)(0) = 1. In Fig.ll, both of the dilaton field and the 
conformal factor a slowly increase monotonically but the scalar curvature 
vibrates rather quickly like trigonometric functions but the amplitude slowly 
increases. In Fig. 12, where we chose A^ = 100 and the initial condition 
0(0) = ct(0) = a(3)(0) = 0, 0(0) = (t(0) = o-(O) = 1, the conformal factor p 
behaves in a rather complicated way, i.e., p increases at first but decreases 
after that with slow and small oscillations. On the other hand, increases 
monotonically and reaches a singularity in finite conformal time. Since the 
conformal factor p does not become positive infinity, the finite conformal time 
corresponds to finite cosmological time. The 4d scalar curvature vibrates and 
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the amplitude increases quickly near the singularity. The 4d scalar curvature 
oscillates slowly in Fig. 13, where = 100 and the initial condition is 0(0) = 
cr(0) = (t(0) = a(0) = a^^\0) = 0, 0(0) = 1. The conformal factor p increases 
at first but decreases after that with very slow oscillations and increases 
monotonically and reaches a singularity in finite conformal and cosmological 
times. In Fig. 14, where N = 100 with the initial condition 0(0) = 0(0) = 
cr(0) = d{0) = a^^\0) = 0, (j(0) = 1, the scalar curvature oscillates very 
quickly and the period becomes short with the passing of conformal time. In 
this solution, the dilaton field also oscillates but the amplitude becomes 
small and the conformal factor p slowly and monotonically increase. 

Hence, unlike the case of pure Einstein gravity with conformal matter 
back-reaction on gravitational background (see ref.[^]) it is more difficult to 
realise the inflationary Universe in the theory with dilaton under considera- 
tion. 



7 Discussion 



In the present work we studied quantum cosmology for some 2d and 
4d dilatonic SGs with dilaton coupled matter. Working in the large ap- 
proximation we demonstrated that quantum matter back reaction does not 
really help to solve the singularity problem unlike the case with no external 
dilaton. Of course, this does not mean that we claim that the singularity 
problem cannot be solved as result of quantum fluctuations. Rather our re- 
sults indicate that the solution of this problem perhaps could be possible only 
in complete quantum SG models, i.e. taking into account quantum effects 
of all fields (graviton,dilaton,gravitino,...). It would be really interesting to 
study such a problem but it is not easy to do since we do not know the 
correspondent expression for effective action (at best, it is known only as an 
expansion in the curvature, see [|1^ for a review). Note, nevertheless, that we 
found few non-singular Universes which are non-singular at early times but 
become singular at late times (or vice- versa). 

Another interesting aspect of our work is the following one. If we exchange 
the role of time and radius for many of our numerical solutions we find a static 
object which may be interpreted as a 2d black hole (BH). A few years ago 
there was a lot of activity in the study of quantum dilatonic 2d BHs (for 
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very incomplete list of references see and refs. therein). Hence, using 



our approach one can study 2d dilatonic BHs for the models of dilaton SG 
with quantum matter. It is also interesting to note that recently discussed 
effect of BH anti-evaporation p4[ also takes place in above BHs realized in 



SGs as it could be confirmed by easy calculation. (It is natural as only the 
coefficient of the Polyakov term in the anomaly induced effective action is 
changed due to the fermion's contribution). 

As one more possible generalization of this work one can consider not 
only conformally flat but more complicated FRW cosmologies for 4d theory. 
However, we do not expect that the results will be qualitatively very different 
at least in the models of the same sort. The reason is that main effects come 
through time dependence which will be similar. 
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Figure Captions 

Fig.l (solid line), p (broken line) and i?/100 (dotted line) for N ^l,A^ 
and = 1. 

Fig.2 0/10 (solid line), p (broken line) and R/im (dotted line) for = 10, 

^2 = and A2 = 1. 

Fig. 3 (p (solid line), p (broken line) and i?/30 (dotted line) for iV = 10, = 1 
and = 0. 

Fig.4 (solid line), p (broken line) and i?/500 (dotted line) for N = 100, 

= 1 and A2 = 1. 

Fig. 5 (/) (solid line), p (broken line) and i?4/100 (dotted line) for N = 1, 
A^ = and A^ = 1. 

Fig. 6 (f) (solid line), p (broken line) and i?4/100 (dotted line) for N — 10, 

A^ = 1 and A2 = 0. 

Fig. 7 (f) (solid line), px 10 (broken line) and 9 (dotted line) for N — 1, — 
and = 1. 

Fig. 8 (p (solid line), p x 4 (broken line) and 6 (dotted line) for N = 1, A'^ = 1 
and = 0. 

Fig.9 0x5 (solid line), p (broken line) and 9 (dotted line) for A?" = 10, ^4^ = 1 

and A^ = 1. 

Fig.lO 0x 100 (solid line), pxlO (broken line) and 9 (dotted line) for A^ = 100, 
A^ = 1 and A^ = 0. 

Fig. 11 (solid line), a (broken line) and 4d curvature divided by 100 (dotted 
line) for AT = 1 with the initial condition 0(0) = a{0) = 0, 0(0) = &{0) = 
<t(0) = a(3)(0) = 1. 

Fig. 12 (solid line), o" x 5 (broken line) and 4d curvature divided by 10 
(dotted line) for A^ = 100 with the initial condition 0(0) = cr(0) = a^^^O) = 
0, 0(0) = ct(0) = (t(0) = 1. 

Fig. 13 (solid hue), cr x 10 (broken line) and 4d curvature divided by 10 
(dotted hue) for A^ = 100 with the initial condition 0(0) = a{0) = &{0) = 
a(0) = a(3)(0) = 0, 0(0) = 1. 

Fig. 14 X 500 (solid line), a (broken line) and 4d curvature divided by 5 
(dotted hue) for N = 100 with the initial condition 0(0) = 0(0) = a{0) = 
a{0) = a(^\0) = 0, (7(0) = 1. 
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